Escape time statistics for mushroom billiards 
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Chaotic orbits of mushroom billiards display intermittent behaviors. We investigate statistical 
properties of this system by constructing an infinite partition on the chaotic part of a Poincare 
surface which illustrates details of chaotic dynamics. Each piece of the infinite partition has an 
unique escape time from the half disk region, and from this result it is shown that, for fixed values 
of the system parameters, the escape time distribution obeys power law l/t^ c . 

PACS numbers: 05.45.-a 



I. INTRODUCTION 

Fully chaotic dynamical systems such as the baker 
transformation and the Arnold's cat map are statisti- 
cally characterized by, for example, exponential decay 
of correlation functions with decay rates given by the 
Pollicotte-Ruelle resonances (See Ref. Q and references 
therein) and exponentially fast escape from regions of 
phase spaces with the escape rate given by the posi- 
tive Lyapunov exponents and KS (Kolmogolov-Sinai) en- 
tropy [J y[ . These properties are outcomes of the uni- 
form hyperbolicity, which means the uniform instability 
of chaotic trajectories. 

In contrast to such ideally chaotic systems, phase 
spaces of generic Hamiltonian systems consist not only 
of non-integrable chaotic regions but also of integrable 
regions (torus), where motions are quasi-periodic and 
therefore the uniform instability may be broken in these 
systems. In fact, generic Hamiltonian systems frequently 
exhibit power law type behaviors, which is due to oc- 
casional trappings of chaotic orbits in neighborhoods of 
torus. Although these phenomena are observed in many 
systems 0, H S B 0|, analytical derivations of decay- 
ing properties of correlation functions and sticking time 
distributions are difficult mainly because there exist com- 
plex fractal torus structures. 

In order to understand power law behaviors in dy- 
namical systems, non-hyperbolic 1-dimensional map- 
pings have been studied by several authors (e.g., Refs. 
US El E El El)- Therefore it is natural to imag- 
ine connections of these non-hyperbolic maps and mixed 
type Hamiltonian systems, however, extensions of these 
maps to 2-dimensional area-preserving systems are un- 
known (but see Refs. [ID, [HI). Thus it is important 
to elucidate the properties of non-hyperbolicity which is 
typical in the mixed type Hamiltonian systems. 

The mushroom billiard, which has been proposed by 
Bunimovich recently [TtJ , is expected to be a candidate of 
analytically tractable model for such problems of mixed 
type systems. This is because the mushroom billiard sys- 
tem docs not have the fractal torus structures and chaotic 
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FIG. 1: The shape of the board of the mushroom billiard 
(the solid lines), which consists of a half disk (the hat) and a 
rectangle (the foot). A point particle inside the board moves 
freely except for the elastic collisions with the walls. The 
absolute value of the angular momentum \L\ equals to the 
distance between the origin and the trajectory (See Sec. II A.) . 
The boundaries L (9) of the first escape domain L~ (8) < L < 
L + {6) are also displayed (See Sec. II B.). 

and torus regions are sharply divided. (See Ref. [II], for 
other example of such systems.) Thus the mushroom 
billiard system can be thought as an ideal model for un- 
derstanding mixed type Hamiltonian systems, and it has 
already been under active researches [IH, [2(| [2l|, IHJ . 

In this paper, we give a theoretical derivation of the 
escape time distribution for fixed values of the system 
parameters. In Ref. [l9j | . it has already been shown nu- 
merically that it obeys power law, and our result agrees 
theirs perfectly. In order to derive the escape time dis- 
tribution, we begin with the construction of an infinite 
partition on a Poincare surface, which reveals detailed 
dynamics in neighborhoods of the outermost tori. 

This paper is organized as follows. In Sec. II, we intro- 
duce the mushroom billiard system, and define a Poincare 
map and its inverse transformation. In Sec. Ill, we con- 
struct the infinite partition by using the inverse of the 
Poincare map recursively. And in Sec. IV, the escape 
time distribution is derived from the structure of the in- 
finite partition. A brief discussion is given in Sec. V. 
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II. POINCARE MAP AND ITS INVERSE 

The mushroom billiard is defined by the motion of a 
point particle on the billiard board depicted in Fig. [TJ 
This board consists of a half disk (the hat) of radius R 
and a rectangle (the foot) of width r and height h [i~7l |. 
We use the polar coordinate (u, 9) , and the Cartesian 
coordinate (x,y); we set the origin as the center of the 
half disk in both cases. The angle variable 9 is denned as 
the angle between the position vector of the point particle 
and the vertical line (See Fig. [TJ. 



A. The definition of the Poincare Map 



We define a Poincare surface at the arc of the semicircle 
x 2 + y 2 = R 2 (y > 0) with negative momentum of the 
radial direction, namely, just after the collision with the 
arc. For the coordinate of the Poincare map, we use the 
angle 9 and the associated angular momentum L. This 
Poincare map $(£, 9) is area-preserving; it can be proved 
through a direct calculation of the Jacobian of the map 
which equals to I everywhere. This coordinate system is 
slightly different from the Birkhoff coordinate, because 
the former is defined only on the arc, but the latter on 
the whole boundary of the billiard board. 



We also set the kinetic energy as v 2 



I. Although 



this setting is not essential, it is convenient to calculate 
the angular momentum L; the absolute value of the an- 
gular momentum \L\ equals to the distance from the ori- 
gin to the trajectory. In Fig. [TJ for example, if a point 
particle moves on the line PA in the direction described 
in the figure, its angular momentum L (L > 0) equals 
the length of the segment AO (the dashed line), and if 
a point particle moves on the line PB in the direction 
described in the figure, the absolute value of the angular 
momentum — L (L < 0) equals the length of the segment 
BO (the long-dashed line). 

We display an example of the Poincare surface in 
Fig. [21 The Poincare map Q(L, 9) is defined on T> = 
{(L,9) e [-R, R] x [-7r/2,7r/2]}; the region \L\ < r 
is chaotic, and \L\ > r is torus. The Poincare map is 
symmetric with respect to the origin (L, 9) = (0,0), i.e., 
$(L, 9) = — $( — L, — 9). In the subsequent subsections, 
we will restrict the domain of the Poincare map to the 
region of the negative momentum in order to simplify the 
analysis. 



B. The first escape and re-injection domains 

Let us consider a point (L, 9) on the Poincare surface 
such that the original orbit of the billiard system (the 
continuous time flow) starting from this point escapes 
from the hat region to the foot without no collision. We 
define the first escape domain T>\ as all such points on the 
Poincare surface. The boundary of T>i can be calculated 
analytically as follows; fix the angle 9 on the Poincare 
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FIG. 2: The Poincare surface for R = 1, r = 0.5 and h — 1. 
The region \L\ < 0.5 is chaotic and the other integrable. 

surface, if the angular momentum L satisfies the relation 
L-{9) < L < L+(6), then (L,9) G V x , where L±(0) are 
displayed in Fig. [TJ More precisely, (9) are defined as 



i ± (6') = ±r sin(7r — ip±) 

= ±rRcos9 (R 2 +r 2 T^Rrsint 



-1/2 



(1) 



where the angles ip± are defined as in Fig. [TJ and given by 
tany>± = Rcos9/(— Rsin9±r). Next, let us consider the 
the domain with negative angular momentum, L~(6) < 
L < 0; the positive domain L > can be treated in the 
same way because of the symmetry. Solving the equation 
(L~(9)) 2 > L 2 in terms of 6, the first escape domain for 
L < can be represented as 



Pi = {(L,9) e V- | 6-{L) < 9 < 8+{L)}, 



(2) 



where V~ = {(L, 9) e [—r, 0] x [-7r/2, 7r, 2]} is the chaotic 
region with negative angular momentum and ± (i) are 
defined as 



(L) = arcsin < 



-L 2 ± ^L 4 -L 2 {R 2 



r 2 R 2 



rR 



.(3) 



The functions 9 ± (L) defines the boundary of the first 
escape domain T>\. 

An orbit of the billiard flows starting from the first es- 
cape domain exits the hat region and stays in the foot for 
some times; and then it returns to the hat and reaches 
again to the Poincare surface. We define the re-injection 
domain T> m on the Poincare surface as all such just re- 
turning points , more precisely, we define V m :— $(2?i). 
2?i n can be derived in the same way as T>± ; 

V in = {(L, 9) e V- I er(L) < 9 < 6+(L)}, (4) 

where the boundary 9 in (L) is defined by 



L 2 ± ^/L 4 - L 2 {R 2 + r 2 ) 



r 2 R 2 



rR 



.(5) 



In FigH 9 ± {L) and 9f n (L) are displayed for R 
r = 0.5 and h > 0. 
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FIG. 3: The boundaries 0+(L) (the solid line), 0~(L) (the 
dotted line), 6^ n (L) (the dashed line), 9r(L) (the dotted- 
and-dashed line), and 8o(L) (the dotted-and-double-dashed 
line) are described. We set the system parameters as R = 1, 
r = 0.5 and h > 0. The only negative part of the angular 
momentum L < is displayed because the Poincare surface 
is symmetric. 



C. The inverse of the Poincare map 

When an orbit of the billiard flows collides with the 
boundary 

W := { (x, y) | x e [-R, -r] U [r, R], y = 0}, (6) 

the angular momentum changes its sign. We should take 
into account the collisions with this boundary W, be- 
cause we reduce the Poincare map to the domain of neg- 
ative momentum L < 0. Let us consider the domain 
D ~\£> in and its inverse image $~ 1 (D~\T>i n ). The orig- 
inal orbit I (namely flow) connecting a point (L,6) G 
V-\V in and ^(L, 9) E ^(l?" \X> in ) are classified into 
two classes for fixed 9 (See Fig. Ufa).): when L < L n (9), 
the orbit I is a line segment, namely there is no collision 
with the boundary W and when L > Lq(6), the orbit I 
consists of two line segments, namely there is a collision 
with the boundary W. In Fig. HJa), an orbit for the crit- 
ical case L = Lq(6) is displayed by the dashed line. We 
can define L (9) by 



L (9) = -i?cos(^ - 9), 



(7) 



where "0 is defined as depicted in Fig.[?Ja). Furthermore, 
using tan?/> = cos 9/(1 — sinfl), we have 



L {9) = -R(2-2sin9)- 1/2 cos9. 



(8) 



Solving the inequalities L < Lq(9) and L > Lq(9) in 
terms of 9, we have the result: when 9 < 9q(L) there is 
no collision with the boundary W, and when 9 > 9q(L) 



FIG. 4: (a) A classification of orbits for fixed : if L < 
Lo(6) (< 0) there is a collision with the wall W and if Lo(6) < 
L (< 0) there is no collision. Note that we consider only 
L < 0. In the figure, we display the critical case L = Lo (9) by 
the dashed line, (b) The inverse image of the Poincare map, 
which is derived from 8 — 9' + tv — 2<j> and </> — arcsin \L/R\. 



there is a collision with the boundary W, where 9q(L) is 
defined by 



'21? 

%(L) = arcsin ( — - 1 



(9) 



Using these results and definitions, we can construct 
the inverse of the Poincare map <i> _1 on T)~\D m as fol- 
lows (See Fig.Hb)), 



tt — 2 arcsin 



L, 9 — 2 arcsin 



if 9 < 9 Q (L) (10) 
, if 6 > 9 a (L) 



Notice that the angular momentum is unchanged by the 
collisions with the arc, and that we restrict the domain 
of the inverse map on the region L < by identifying 
the points (-L, —9) with (L, 9) when the point particle 
collides with the wall W, i.e., when 9 > 9q(L). 



III. THE INFINITE PARTITION 

Using this inverse map < &~ 1 , we can define the n-th 
escape domain T> n recursively: 



P n = $- 1 (D n _ 1 \O in ), (n = 2,3,- 



(11) 



Note that we should remove T> m from T> n -i in the re- 
cursion relation Eq. (IIip . because the inverse image of 
the injection domain $ _1 (Pi n ) equals to the first escape 
domain T>\ . 
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FIG. 5: (a) The first to fourth escape domains T>i—T>4 are displayed. The second escape domain T>2 is separated two parts. 
The remaining part consists of three pieces, which we define as Eo,E\ and i?2- (b) The regions T\ are displayed. The dotted 
line indicates the boundary of the re-injection domain. 

We restrict the parameters as R = 1, r = 0.5 and h > 
in the following and explicitly derive the boundaries of 
the n-th escape domain T> n . First, we derive explicitly 
the first to fourth escape domains, in order to confirm 
that these four domains fill the domain T>~ like Fig. [5](a) 
except for the three regions E\ , E 2 and E$ . Then the 
boundaries of the n-th escape domain (n > 5) can be 
derived recursively. 

Let us start with n = 2 of Eq. (p~T|) ■ The domain 
Di/T>- m can be divided into three pieces: 

Pi/An = {(L,9) I 9-{L) < 9 < 9o(L)} 

U{(£, 6) | 9 (L) <9< d+(L), -r<L< L~(G)} (12) 

U{0M) I O (L) < 9 < 9r(L), L~(0) < L < 0}. 

where L~(0) is defined by Eq. ((T|) (See also Fig. 03(b)). 
In Eq. (Ti"2"|) , we abbreviate the expression (L, 9) G T>~ to 
the one (L, 9) for simplicity; we use the same abbreviation 
in what follows. Let us represent the three sets on the 



right hand side as T\, T\, and T\ , respectively; namely, 
D i/2?i„ = T\ U T\ U T\ . These three sets are displayed 
in Fig. E{b). Using these notations and the inverse map 
<i> _1 [Eq. <Hnj)], the second escape domain T> 2 is given by 

V 2 = $-!(!?! /V in ) 

(13) 

= r'^jur'^ur 1 ^ 1 ). 

Let us denote the 9 component of <& _1 (L,#) as ^ L (9). 
and inverse image of T\ as T\ = < I> _1 (J r 1 1 ). Using these 
definitions, we have 

T\ = {(L,9) | 9 L (6-(L)) < 9 < * L (0o(L))} 

(14) 

= {(L,0) | 6-{L) +7T- 2arcsin|L| < 9 < §)}, 

where the upper bound for 9 is 9 = tt/2. Note that 
the angular momentum L is unchanged under the inverse 
map <I> -1 . Similarly, we define the inverse image of T\ 
as T 2 = ^^(J 7 ^), and we have 



jq = {(L,9)\* L (9 (L))<9<t> L (9+(L)), -r<L<L-(0)} 

(15) 

= {(£,6>)| - f < 9 < 9+{L) - 2arcsin|L| , - r < L < L~(0)}, 
where the lower bound for 9 is 9 = —tt/2. Finally, let us define the inverse image of J-\ as T\ = §~ X {T\ ), and we get 

jq = {(L,9) I * L (0 O (L)) < 9 < ^ L (0^(L)), L"(0) < L < 0} 

(16) 

= {{L,9) | - f < 6> < 0-(L), L~(0) < L < 0} 

where we have used the relation 9~{L) — = — 2arcsin|i|. In the Eq. (TIB)) , the lower bound for 9 is 9 — —tt/2; 

the upper is 9 = 6~(L), which is equivalent to the lower bound of the domain T>\, These three sets {J-?, J 7 !,^ 7 !} are 
displayed in Fig. [D(b). 
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FIG. 6: (a) The infinite partition constructed in terms of the escape time. The solid lines represents boundaries between regions 
of different escape times. The right hand side of the broken line is the injection domain, (b) A magnification of Fig. H|a) in 
a neighborhood of the outermost tori (L = —0.5). The boundaries for the domains with the escape times longer than 32 are 
omitted. 



Next, we derive the third escape domain T>^. It can be proved that T>2 H T> m = 0, because the relation 9 (L) + 
7r — 2 arcsin |L| > 9^ n {L) holds. (See the lower bound of the domain T\ which is defined by the second line of the 
Eq. Crj}). It follows that X> 3 = $ _1 (X> 2 ) = ^{^D U U by Eq. ([11]), where the three sets of the 

right hand side can be calculated, respectively, as 



S- 1 ^) = {(L,6) 
= {(L,6) 

*-\J%) = {(L,6) 
= {(L,6) 

$"H^!) = {(L,6) 
= {(L,9) 



^ L (9~ (L) + 7r — 2 arcsin \L\) < 9 < * L (f )} 
6>~(£) + 7r-4arcsin|L| < 9 < § - 2 arcsin \L\}, 

* L (-f) < 9 < ^ L {9+{L) -2 arcsin |£|), - r < L < L-{0)} 

§ - 2 arcsin |L| < 9 < 9+{L) + tt - 4 arcsin \L\ , -r<i<i~(0)}, 

* L (-f ) < 9 < ^ L {9-{L)) 1 L-(0) < L < 0} 

f - 2 arcsin |L| < 9 < 9~ (L) + n ~ 2 arcsin \L\ , L"(0) < L < 0}, 



(17) 



where the relation $ H-^i) C An holds, because of the inequalities 9f n {L) > § — 2 arcsin |L| and 9 (L) + tt 



4 arcsin \L\ > 9 in (L). Therefor, the 4-th escape domain V4 is represented by D4 — $ 1 (2?3/2?i n ) : 
where we define 

T\ = {(L,6) I ff+(L) < 6 < 9+ (L)+tt- 4 arcsin \L\, -r<L<L-(Q)} 
jq = {(£,0) I 0+(£) < < 6-{L) + tt - 2 arcsin L~(0) <L < 0} 



*- 1 (^?)u*- 1 (^'!), 



(18) 



These sets {-^1,^2} are displayed in Fig. H^b). Thus, the fourth escape domain V4 is given by the union of the 
following sets: 



$ _1 (^f) = {(L,6) I 9+(L) < 9< 6»+(i) +7T- 6 arcsin - r < L < L~(0)} 
$ _1 (.Ff) = {(L,6l) I 6»+(i) < 6» < 6»-(L) +7r-4arcsin| J L| , £-(0) < L < 0}, 



(19) 



where we have used the relation 9 + (L)~9^ n (L) — —2 arcsin \L\. The lower bounds for 9 of these two sets are 9 = 9 + (L); 
and the upper bounds of the set Q~ 1 (J : 2 ) is equivalent to the lower bound of T>^. 



From the above results, the domain V are covered by the four sets 2?i,2?2,2?3, and V4 except for the three 
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FIG. 7: The cumulative distributions of the escape times for 
discrete time (the solid line) and for continuous time (the 
dotted line) in log-log form. The dashed line is the analytical 
result F(t) ~ 1/t 2 . The inset is a magnification of the discrete 
time case, which shows a clear stepwise structure. 



FIG. 8: Auto correlation function of the absolute value of 
the angular momentum \L\ — (\L\) (the solid line) in log-log 
form, where (■) means the ensemble average in terms of the 
Lebesgue measure in the chaotic domain. The broken line 
represents a function f(n) ~ 1/n, which is a guide to see. 



regions Eq,Ei, and E 2 as illustrated in Fig. OJa) (b) . The 
boundaries of the n-th escape domain (n > 5) can be 
obtained by recursively calculating the inverse mapping 
<i>~ 1 of the upper bound of the set $ _1 (J r j 5 ), which is 
given by Eq. (JT5J) as 9 = 9+(L) + n - 6 arcsin \L\. Thus, 
let us define the boundary between the domains T> 3n+ i 
and £> 3 ( n+1 ) +1 as 6 3n+1 (L) (n > 1), we have 

6» 3n+ i(L) = (9+ (L) + nn-6n arcsin \L\ (20) 

Similarly, defining the boundaries between T>^ n+ 2 and 
^3(«+i)+2 as 9 an+2 (L), and between V 3n and V 3 ( n+1) as 
3n (L) , we have 

3 n+2{L) = 6+{L) +7MT- 2(3n+ I) arcsin |L| , 

(21) 

3n (£) = 0+{L) +nn- 2(3n- 1) arcsin \L\ 

In Figs. [Hfa) and (b), we depict these boundaries up to 
n = 31. 



IV. THE ESCAPE TIME DISTRIBUTION 

Finally, we derive a scaling property of the escape time 
distribution f eS c(n) approximately. The escape time is 
defined by the number of collisions with the arc of the 
semicircle just after an orbit enters the hat region until 
it escapes from there. 

Since the Poincare map is area-preserving, which is the 
universal property of the Poincare map of the Hamilto- 
nian systems, the physically natural invariant measure 
of the Poincare map is the Lebesgue measure. Thus the 
probability that the escape time equals n is given by the 
Lebesgue measure of the region T> n n An, that is 



/ esc (n) = S*. Note that = = (n = 2, 3, • • •) 

and thus only S 3 ^ l+1 (n = 1, 2, • • •) has finite values (See 
Fig. Ha) (b).). 

Let us derive the cross points of the lines 6 3n+ i(L) = 
9 + {L)+n-K-6n arcsin \L\ [Eq. and 6(L) = 0. Using 
the Taylor expansions arcsinx w ±-| + ^(it|) (as 

x w ±i), we have 

it 2 r- 1 ( 22 ) 

-- + ^(3n + l)(2L + 1) + V2(2L + 
o v 3 

asm oo. Setting this to 0, and solving in terms of L, 

we have 2L + 1 ~ „ 1 [ , ~ —. Thus we find the width of 

3n+l n 

the (3n + l)-th escape domain is proportional to 1/n 2 . 
It follows that the area of the n-th escape domain S n 
behaves as 

S n ~ \, (23) 

as ti — ► oo. Note that S 3n +i = S 3n = S 3n -i for 
n = 2, 3, • • •, because the domains 2?3„ and I?3„_i has no 
intersection with the domain V ln . The Eq. (f2"3"]l means 
the fact that the partition constructed in the previous 
section is infinite. Finally, we get 

S 3n +l = S 3n+ \ — 53„ + 4 ~ — 3 , (24) 

as n -» oo. And, as mentioned above, the equations 
S 3 n n = S 1 3 r n _ 1 = hold. This power law perfectly agrees 
with the numerical results shown in Fig. where the 
cumulative distribution 

oo 

i^ sc (n):= ( 25 ) 

j=n+l 
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is plotted by the solid line. Note that this numeri- 
cal result have already been reported by Altmann et. 
al. 19]. In the inset, a magnification is displayed, which 
shows a clear stepwise structure with decreases exactly 
at n = 3k + 1 (k — 1, 2, • • •). This implies that / C sc(3fc) — 
/csc(3/c + 2) = 0, and that the only / esc (3fc + l) have finite 
values (k > 2). Thus, these results also agree with the 
analytical results. 

V. CONCLUDING REMARKS 

In conclusion, we have derived the escape time distri- 
bution by constructing the infinite partition in terms of 
the escape times. Note that, however, the escape 'time' 
in this paper is the number of collision until the particle 
escapes. Thus the escape time of the continuous time 
flow might be slightly different from ours. But the scal- 
ing exponents should be the same, because the flight time 
of the chaotic orbits between collisions in the hat region 
is non-vanishing. This is confirmed numerically and the 
result is displayed in Fig. [7] which shows the agreement 



of the scaling exponents of these two distributions. 

There are several points that should be verified in fu- 
ture studies. First, the correlation functions of this sys- 
tem exhibits power law behavior (Fig. El see also [17j). 
and it is expected that there are relations between the 
scaling exponent of the escape time distributions and the 
correlation functions. Second, it is important to elucidate 
whether the results in this paper is general or not for 
other parameter values R ^ 2r. Third, it is also impor- 
tant to consider which property of the mushroom billiard 
is the universal features of the mixed type Hamiltonian 
systems. 
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